On the Ordering of Na + Ions in Na x Co02 
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Abstract. The influence of electrostatic interactions on the ordering of sodium ions in Na A Co02 is studied theoretically 
through Monte-Carlo simulations. For large x small di- or tri-vacancy clusters are stable with respect to isolated Na vacancies. 
At commensurate fillings these small clusters order in triangular superstructures. These results agree with recent electron 
diffraction data at x = 1/2 and 3 /4. We have performed neutron Laiie diffraction experiments at higher x, which confirm the 
predictions of this simple model. The consequences on the properties of the electronic charges in the Co layers are discussed. 



In the past decade, Na A Co02 has emerged as a sys- 
tem of fundamental interest because of its high ther- 
moelectric power, unusual magnetic properties and pos- 
sibility, when hydrated, to superconduct 1 1 ] . The den- 
sity x of intercalated sodium atoms tunes directly the 
density of states of quasi-particles in the metallic CoO? 
sheets. For simple fractional fillings, electron diffraction 
(ED) experiments 13| provide evidence for the order- 
ing of Na + ions at room temperature. This non-trivial 
order is generally believed to be the consequence of in- 
tricate interactions between Na + ions, mobile electrons, 
and phonons. We show here that the dominant process 
driving this ordering is simply the Coulomb potential be- 
tween Na + ions. Sodium ions can occupy sites between 
oxygen atoms, which form two interpenetrating triangu- 
lar lattices denoted Nai and Na2. Cobalt ions lie above 
and below Nai sites, resulting in an extra energy cost 
A = A c + A sr relative to an Na2 site (A c in the long-range 
Coulomb part and A sr a short range repulsion between 
Na and Co ions shells). Na ions are bigger than the dis- 
tance between nearest Nai and Na2 sites, forbidding si- 
multaneous occupancy of nearest-neighboring sites and 
that makes Coulomb energy minimization highly non- 
trivial. Our model Hamiltonian includes: (i) a long-range 
Coulomb potential, e/(4neeor) (the dielectric constant, 
taken as isotropic, is fixed at e = 6 to account quanti- 
tatively for the variations of the chemical potential, in 
terms of x, measured in cell devices |4]); (ii) Na-Na ion- 
shell repulsion V=0.04e V for neighbors on the same sub- 



lattice; and (iii) on-site energy A S7 - = O.OleV. A system 
of two sheets of Na ions, each containing a maximum 
number of 1176 ions, intercalated between two C0O2 
layers (with periodic boundary conditions in 3 dimen- 
sions) is simulated at finite temperature and fixed chemi- 
cal potential fj, through a Grand-Canonical Monte-Carlo 
method. The hexagonal cell parameters a=2.84A and 
c=10.87A are fixed to those measured in the x = 0.75 
phase, and elastic deformations are neglected. We main- 
tain neutrality at each step by varying the charges in the 
Co layers assuming a uniform spreading. 

The organization of Na ions is driven by the sponta- 
neous formation of multivacancy clusters, as illustrated 
in Figure 1 . At large distance d, the repulsion energy be- 
tween two vacancies decreases as 1 jd. However neigh- 
boring vacancies can reduce their energy by promotion of 
a Na2 sodium to the central Nai site. The resultant clus- 
ter has a net charge 2e~ spread over three sites and sub- 
stantially lower energy, as the central sodium is now fur- 
ther from its neighbors. Formation of a tri-vacancy clus- 
ter follows a similar process and three vacancies combine 
with three Na2 — > Nai promotions. 

Monte-Carlo results at room temperature are presented 
in Figure 2. The Na concentration versus chemical poten- 
tial is characterized by plateaus at x — |, |, |, |, com- 
parable to those appearing in electrochemical cells 0]. 
The the j -plateau corresponds to an {a'=2a-b, b'=2b} 
ordered structure with equal occupation of Nai and Na2 
sites in perfect agreement with ED spectra |2]. It can be 



described as a dense lattice of divacancy clusters. The I- 
structure is a {a'=2a+b, b'=-a+3b} triangular array of 
divacancies. At x = f a triangular array {a' = 2a + 2b, 
b'= -2a+4b} of trivacancies is obtained, in agreement 
with ED |3|3[. No plateau is observed at higher x, but 
short-range ordered random mixing of di- trivacancies 
appears in the simulations. Nevertheless, at x > 0.8 many 
commensurate triangular arrays of di- trivacancies are 
possible and we have therefore investigated experimen- 
tally this range though neutron diffraction. 

Single crystal data were collected on the SXD diffrac- 
tometer at the ISIS pulsed neutron source, Rutherford 
Appleton Laboratory (UK) with area detectors covering 
a solid angle of 2n steradians. At x = 0.92 ± 0.03, Fig- 
ure 3 shows superstructure peaks appearing on a circle of 
radius q « 6A6nmT i centered on the main Bragg points. 
These are compatible with a triangular array of divacan- 
cies with unit-cell vectors {a'=3a+b, b'=-a+4b}. 

Sodium ordering induces a periodic Coulomb poten- 
tial in the Co layers. The depth of the electrostatic po- 
tential wells, «100meV, is substantially larger than the 
quasi-particle hopping frequency «10meV and so will 
localize holes. While total energy differences are in- 
significant (wO.lmeV), the Coulomb landscape varies 
drastically in terms of position of vacancy clusters in 
successive layers along c. For x — \ there are only 
two configurations. Figure 4 (a) corresponds to divacan- 
cies in successive layers displaced by b. It has a one- 
dimensional character with high (dark) and low (light) 
Coulomb-potential stripes in agreement with the low 
and high spin charge-ordered stripes found using neu- 
tron diffraction |5]. A large magnetic field parallel to the 
(a,b) plane favors alignment of holes (divacancies) in the 
c-direction corresponding to Figure 4 (b) with 2D char- 
acter and conducting Co ions at intermediate valencies 
(gray circles). This could explain the abrupt change from 
an insulating antiferromagnet to a conducting state with 
reconstruction of a 2D Fermi surface when a field of 40T 
is applied at T=50K |6]. The potentials corresponding to 
the x — | and | structures are discussed in Ref. Q. 
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FIGURE 2. Na concentration x (solid line) and relative pro- 
portion of Na2 sites (dashed line) versus chemical potential jl 



In conclusion, at x > |, the ordering of sodium ions 
governs the electronic and magnetic properties of quasi- 
particles (holes) in the Co layers. Pure electrostatics with 
non-trivial constraints on a bipartite lattice leads to va- 
cancy clustering. For simple fractional fillings di- and 
trivacancy clusters order long range. 
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FIGURE 1. Di- and tri-vacancies. Dark (light) circles corre- 
spond to occupied Na2 (Nai) sites 
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FIGURE 3. Neutron spectrum in (h,k,6) plane at x fa 0.92 



FIGURE 4. Na Coulomb potential in Co planes at x = 5 with 
divacancies aligned along c axis (a) and displaced by b (b). The 
white (black) circles represent Co~ 3+ (Co~ 4+ ) ions. The gray 
circles represent conducting Co ions at intermediate valencies. 



